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Abstract 

In this paper, we give a few results on the local behavior of harmonic 
functions on the Sierpinski triangle - more precisely, of their restriction to 
a side of the triangle. First we present a general formula that gives the 
Holder exponent of such a function in a given point. From this formula, 
we deduce an explicit algorithm to calculate this exponent in any rational 
point, and the fact that the derivative of such a function is always equal 
to 0, oo or undefined. 

1 Introduction 

1.1 Notations and conventions 

In all this text, uj stands for the cubic root of —1 with positive imaginary part. 

For z e C, we define h z : C —¥ C to be the homothety of center z and ratio 
\: =»±». 

N stands for the set of all nonnegative integers. 

When talking about a positive quantity, we shall say that it is positively 
bounded if it is bounded from above by a finite constant and from below by a 
positive constant. This is equivalent to saying that its logarithm is bounded. 

When we say that some derivative is well-defined, we shall mean that the 
rate of change has a limit in R, i. e. we shall also implicitly admit infinitive 
derivatives. 

We will often need to make estimations up to a positively bounded multi- 
plicative constant. As all norms on a finite-dimensional vector space are equiva- 
lent, in such a context, we have no need to distinguish them. So we shall usually 
simply write "|| • ||", implying that the statement we are making is true for any 
norm. However, to perform some calculations, we will need to use some specific 
norms: we shall then use an index to specify which norm we are talking about. 

We will also often need to manipulate products of the form F ai . . . F an , 
where F is some kind of operator that can take one of two values F and F\ 
(for example h, M, M, M and so on). We shall then simply write F ai an , or 
even F w (if we have already agreed that w = a± . . . a n ), as shorthand for such 
products. 

Everything written in positional notation with a radix point is assumed to 
be in binary (for example \ =0.1), unless the context clearly shows otherwise 
(for example 7r ~ 3.14 . . .). 

As everyone knows, dyadic rational numbers have two binary expansions. 
Being given s G [0,1) (resp. s G (0,1]), we call its upper (resp. lower) binary 



expansion the unique binary expansion of s that does not end with an infinite 
sequence of l's (resp. of O's). For the problems that we will study, both ex- 
pansions will play perfectly symmetrical roles; hence we will often simply say 
"binary expansion" without specifying, and we shall mean "upper or lower binary 
expansion". 



1.2 Statement of the problem 

Let u be the unique continuous fonction from [0, 1] to 1Z 3 that satisfies: 

u(0) = , u(l) 

u o hf) — Mq o u 
u o hi = M\ o u, 

% 2 2 N 



(1) 



where we set 



Mn 



l\ 2 1 
\0 1 2, 



Mi 




(see also Figure [5]). We shall justify later the existence and the unicity of this 
function. 

In [lj, Kirillov formulates the question (called "Problem 2", near the end of 
Section 3.2) : 

Compute explicitly the derivative u'(t) whenever it is possible (e.g. 
at all rational points). 

The purpose of our paper is to answer this question. Note however that our 
notations are slightly different from Kirillov's: what he calls u corresponds in 
fact to a projection of the function that we have decided to call u. Thus we 
shall actually deal with a slightly more general question. 



1.3 Plan of the paper 

We start, in section[2] by explaining the interest of the function u. While mostly 
following Kirillov's book, we construct harmonic functions on the Sierpinski tri- 
angle, and we show that u appears as the restriction of some kind of "universal" 
harmonic function to a side of the triangle. The rest of the article is logically 
independent from this section, except for a few notations and definitions. (For 
a more thorough introduction to harmonic functions on the Sierpinski triangle 
and on other fractals, see also Strichartz's book [2].) 

In section [3] we establish a few general results about that function. 

In the following, we study the behaviour of u(t) — u(s) as t — ¥ s. At first 
(section [1]) , we completely describe the behaviour of the direction of this vector. 
This gives the general shape of the curve u(t), and allows us, a few pages further, 
to link the local properties of the functions that Kirillov studies with those of 
u. 



2 



In section [5] we study the norm of this vector. The key result of this paper 
is Proposition [9] that links the asymptotic behaviour of this norm with that of 
an infinite product of matrices indexed on the binary expansion of s. 

In section |5.2| we answer the second part of the question: we present an 
algorithm that allows to calculate the derivative of u (hence also that of u a c b ) 
in every rational point. 

In section |5.3| we answer the question in the general case: we show that the 
derivative of u. when it exists, can only take the values or oo, and that it is 
almost surely equal to 0. 

Finally, in section [5T4| we give some numerical results, which help us make 
a few estimations. We establish a simple sufficient condition for the derivative 
to be equal to 0, and we conjecture one for the derivative to be equal to oo. 
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2 Preliminaries : harmonic function on the Sier- 
pinski triangle 

Definition 1. Let us first define successive approximations of the Sierpinski 
triangle: we start with 

So := {0,1, u}, 
and we define recursively, for all n £ N, 

S„+i := h (S n ) U hi(S n ) U h u (S n ). 

Let us give a name to the union of all these approximations: 

Soo := I J S n , 

and the Sierpinski triangle is then defined as the closure of this union: 



Note than we can naturally see S„ as a graph: So is a complete graph on 
three vertices, and as edges of S n +i, we take the images of the edges of S n by 
the three homotheties. Then it is easy to see that every vertex except 0, 1 and lo 
has exactly four neighbours. (See Figure [I]) In the following, E is a real vector 
space, and we will write "s ~ t" for "s is a neighbour of i". 

Definition 2. Let n £ N. A harmonic function on S„ is a function / : S n — > E 
such that: 

v s eS„\{o,i )W }, = £ £/(*). 

We impose no condition on the values /(0), /(l) and f(u): we consider them 
to be boundary conditions. 
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Figure 1: Discrete approximations of Sierpinski triangle. 



Definition 3. A harmonic function on S is a continuous function / : S — > E 
whose restriction to every S n is harmonic. 

Proposition-definition 1. Let a, b 1 c £ E. Then there exists a unique function 
f that is harmonic on S and that satisfies the boundary conditions: 

/(0) = a, /(l) = b, /(«) = c. 

(This means in particular that the set of all real-valued harmonic functions on 
S is a three- dimensionnal vector space.) We will call this function f a c b . 

Proof. To check this, we will explicitly construct, by induction, the family of 
functions that verify these conditions on successive approximations of iS. 

First we define the (obviously unique) harmonic function on So = {0, 1, tu} 
that satisfies the boundary conditions: 

fo : i y d, 1 i y 6, to i y c. 

Now suppose we have a harmonic function /„ : S n — > E satisfying the 
boundary conditions. We will now define a function f n +i that extends /„ to 
<S n +i- Let s, t, u be any three adjacent vertices of S n that form a triangle having 
the same orientation as Sq. Then 5 n +i still contains these vertices, but also the 
vertices ^^p, ^jjp, (see Figure pi. Conversely, it is easy to see that all vertices 
of <S n +i can be accounted for in this way. We then set 

2f n (s) + 2f n (t) + f n (u) 
5 

2f n (s) + f n (t) + 2f n (u) 
5 

fn(s) + 2f n (t) + 2f n (u) 
5 

Being an extension of /„, this function still satisfies the boundary conditions. 
Checking that it is harmonic on is a straightforward calculation; it is 

enough to do it in s and in , since any vertex is similar to one of those two. 

• The neighbours of s in S n+ i are ^jp, and s ^-, where t' and v! 



f n+ i m = 

fn+1 (^) = 
fn+1 = 



4 



Figure 2: Subdivision of S n . 
are the two other neighbours of s in S n (see Figure [2]) . We have: 

= 2f n (s) + 2f n (t) + f n (u) 2j n (s) + f n {t)+2f n {u) 
5 5 
2f n (s) + 2f n (t') + f n (u') , 2/ n («)+/ n (f) + 2/ n (u / ) 



5 5 

8/n(«) + 3(/ n (t) + / n (M)+/ n (tO+/ n K)) 

5 

8 + 12 , , . 

— /n(s) 

4/„+i(s) 



The neighbours of are s. and i. We have: 

fn+l(s) + fn+ltm + /n+l(^) + /„+l(t) = 



2/ n (a) + / n (t) + 2/ n ( U ) / n QQ + 2/ n (t)+2/ n (tt) 
/4 s ) H g 1" g 1" /«(*) 

8/„(«) + 8/„(t)+4/ B (u) 



= 4/„+i(^) 



Now let us show that these functions are unique. Indeed, consider two harmonic 
functions on S n that satisfy the same boundary conditions. This means that 
their difference vanishes at 0, 1 and lj, and is harmonic as well. Then the 
maximum principle guarantees that it is identically equal to zero. 
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This allows us to define f a c b on Soo, and the above calculation gives us 
an explicit recurrence relation to compute it in any point. Let us rewrite this 
relation in a more convenient way: 

2a + b + 2c 

fa b° h-0 = fa 5 2„+2i,+c 
5 

q + 2b+2c 

fa b ofl l = f 2„+2b+c 5 (2) 

/a C b °h u = /2a + t ,+2c C a + 2i) + 2e . 

5 5 

To extend this function continuously to the whole Sierpinski triangle, we need 
to check that it is uniformly continuous. Let 

D{a,b,c) := max (||o — — c||, \\a — c\\). 

By the maximum principle, we have, for all values of a, 6, c and for all s, t € iSoo> 

||/a C 6(s)-/ C 6 (*)|| < r>(o,6,c). 

On the other hand, it is easy to check that we have 

D{a, i«±|±2c, 20+26+C) < |£,( a) {, )C ). 

Using ([2]), by induction, il follows that if T is an "elementary triangle" of level 
n, we have 

VM€T, ||/a C b ( S )-/a C b (t)|| < (ffC 

where C depends only on a, b and c. Now consider two points of whose 
distance is at most (J) ; then it is easy to see that the respective elementary 
triangles of level n — 1 that contain them are either coincident or adjacent. It 
follows that 

vmgsoo, \s-t\ < (if =► n/a c h ( S )-/ Q c b (<)ii < arc, 

where C is still some positive real constant that depends only on a, b and c — 
which leads to the conclusion. This shows that / Q c fc is indeed well-defined on 
S. ' □ 

Very informally, f a c b describes the shape that a Sierpinski triangle made of 
rubber would assume if it were stretched between three nails fixed in a, b and c 
(see Figure [3]). 

Note that by continuity, the relations |2| are still valid in S. They are very 
important: they supply an intrinsic definition of the functions / Q c b . 

From now on, following Kirillov, we restrict ourselves to a side of the Sier- 
pinski triangle: we let u a c b — / a C b|[ 11 ■ We may then drop the third line from 
([2]). Moreover, Kirillov is mainly interested in real- valued functions, in particu- 
lar four of them: tfi = u ° 1 , ip = Uq 1 \, X — u o _1 i an d £ = Uq 2 ± (in his notations, 
u actually stands for some one function among these four). However, they have 
a drawback: it is difficult to study one of this functions independently from 
others, because (pi forces us to change the values of a, b and c. 
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This leads us to introduce the function u eo eiJ ei , where (eo,ei,e w ) is the 
canonical basis of M 3 . Its main interest is that it contains all the information 
about all the functions u a C b- Indeed, for all a, b and c, we have u a c b — $ou eo e " , 
where $ : M 3 — > E is the linear map that send the canonical basis to (a, b, c). 
(If a, b and c are real, it is simply the linear form with matrix (a b cj.) The 
relations ^ then give us a functional equation on u eo e " ei — which is none other 
than 0. Thus we see that the function u that we have defined in the beginning 
is just a shorter notation (that we shall adopt from now on) for u eo e " ei . Its 
existence and unicity follow from the proof of Proposition-Definition [T] 

3 General properties of u 

Note that all that we shall say in this section can very easily be generalized to 
all of S, that is to / eo e " ei - To do this, it is enough to introduce the matrix 
M u , analogous to Mg and Mi, and to see that a generic point of S can be 
described by an infinite sequence of symbols 0, 1 or w, by analogy with the 
binary expansion. 

We will need to diagonalise the Mj. Any of these matrices has eigenvalues 
1, | and |, with respective eigenvectors ej, Uj and Wi, where: 





Note that 

fx o + fo i + /o o = f^i = 1) 

which shows that our picture is in fact only two-dimcnsionnal: all the values of 
/e e " ei i hence also all the values of u, actually lie in the affine plane 

n := {(x,y,z) eR 3 \ x + y + z = 1} . 

Of course, it is stable by the maps Mj, and the latter induce some affine maps 
on the former (see Figure [i}. In fact, this is why we have chosen to write 
instead of e^: the are elements of the affine plane T-L, namely — under this 
interpretation — fixed points of the affine maps Mj. 

Let us now try to write down an explicit formula for u(s). We first need to 
express s in terms of 0, of 1 and of the hi, which may be done by writing: 

s= lim h ai .,. a (0) = lim h ai ,„ a (1), 

n— foo n— foo 

where s = 0,aia2 ... is a (lower or upper) binary expansion of s. Using the 
above relations, we get the following: 
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Figure 4: Eigenvectors of the operators M*; the plane of the sheet of paper 
corresponds to the affine plane H. Thus the origin of M 3 lies beneath the centre 
of the figure. The triangle K is defined a little bit further; it is added here to 
better illustrate the action of the Mj on H. 



Proposition 2. Let K := H n R> 6e the (full) triangle with the canonical 
coordinate vectors as vertices. Then for all s € [0, 1] and for any initial value 
uq € K , we have 

u(s) = Km M ai ... an (uo), (3) 

n— >oo 

where s — Q.a\a2 ... is iis (lower or upper) binary expansion. 
(See Figure[5]) 

Proof. Let us first rewrite this formula in such a way as to make it clear that 
the result does not depend on the initial value. It is easy to check that K is 
stable by both M and Mi: thus the sequence M ai ... a (K) is decreasing. In 
the proof of Proposition-Definition [T] we have already seen that the diameters 
of these sets tend to 0. It is clear that they are closed. Thus their intersection 
contains a single point: 



o 6 K, { lim M 01 ... 0w (uo)} = H M 01 ... 0n (K). 

Ln— >oo J 1 1 



Vtti 



On the other hand, a trivial induction shows that for any finite word w on the 
alphabet {0, 1}, ^(0) = O.w, where O.w is seen as a binary fraction, and that 
uo h w = M w o u. Hence, by the relations ([!]) and by continuity of u, 



lim M a a 

n— too 
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and hence the formula. □ 
Corollary 3. u is injective. 

Proof. Let (ai)i>\ and (a^)i>i be two sequences with values in {0, 1}. Let k be 
the first index where these sequences diverge, i. e. suppose a\ = a' 1; . . . , au—i = 
a! k _ 1 , ak ^ a! k . Without loss of generality, we may actually suppose that ak = 
and a' k — 1. Choose some uq € K. Then the following statements are equivalent 
(remember that the Mi are bijective): 

lim M ai ... an (u ) = lim M a , a , (u ) 

n—>oo n—too L 

oo oo 

n -Wa,...a„(A')= f] M <-<W 
n^O n=0 

oo oo 

fl M 0ak+1 ... an (K) = fl M 0<+i „ M , n (K) 

n—k n—k 



On the last line, each side is a subset of the respective Mi(K). But a quick 
calculation (or a quick look at Figure |5| shows that the two images of the 
triangle intersect at a single point: 



M Q [K)r\M x {K) = 




= {M (e 1 )} = {M 1 (e )}. 
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Figure 6: Main objects of the vector plane H. 



Thus the previous equality holds if and only if 

' oo 

n M ak+1 .. Mn (K) = { ei } 

J n—k 

\ oo 

p| M 4+i ... < (if) = {e }, 

^ n—k 

which is true if and only if for all i > k, ai — 1 and o! i = 0. In other words, both 
sequences yield the same value if and only if we have J2iLi a i^ 1 — J2iLi a 'i^~ l - 
The "only if part proves our claim, and the "if" part offers additional support 
for the equivalence of different binary expansions. □ 

4 Local behavior of u: existence and direction of 
tangent 

We will now study the behaviour of the vector u(t) — u(s) ad t — > s. In this 
section, we will thoroughly describe the behaviour of the norm of this vector: 
quantitatively through the formula Q , and qualitatively in Corollary [6] 

The difference u(t) — u(s) lives in H := {(x, y, z) e M. 3 | x + y + z = 0}, the 
vector plane parallel to the affine plane H. This vector plane is of course stable 
by Mo and Mi; it is thus natural to introduce the restrictions of these operators 
to %. that we shall call Mq and M\. Each M, has eigenvalues | and |, with 
respective eigenvectors Vi and Wi. Note that in both cases, Vi and Wi actually 
form an orthogonal basis of H. 

Let us now introduce the set 

K := {av + bvi | a, b > and (a, b) ^ (0, 0)} 
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(see Figure |6j) This notation should not lead the reader to think that this set 
is somehow derived from K: it is not, at least not directly. But it does play an 
analogous role, and the following lemma will show precisely in what sense: 

Lemma 4. The set K is stable by the Mi. For all s,t € [0, 1] such that s < t, 
u(t) - u(s) e K. 

Proof. The first part can be checked by a straightforward calculation (see also 
Figure |6| . 

As for the second part, let us first prove the result for dyadic rational values 
of s and t (such that < s < t < 1), by induction on their dyadic valuations. 
Let us write /„ := S n n [0, 1] = | < k < 2™}, so that we have I Q = {0, 1} 
and I n +i = ho(I n )Uhi(I n ). If s, t G Iq, then s = 0, t = 1 and the conclusion can 
be checked by a straightforward calculation. Now suppose that the conclusion 
is true for all s, t € /„; let s < t be two points in I n +i- We need to distinguish 
three cases: s<t<^,s<^<t and | < s < t. We shall only treat the 
middle case, which is the hardest one; the other two are analogous. We may 
write s = ho(S) and t = hi(T), so that we have: 

u(t)-u(s) = u{t)-u{\) + u(|)-n(s) 

= u(fci(T))-«(/n(0)) + u(h (l)) - u(h (S)) 
= Mx(u{T) -u(0)) + M (u(l)-u(S)). 



Applying, in order, the induction hypothesis, the fact that K is stable by the Mi, 
and the fact that K is stable under linear combination with positive coefficients 
(which follows immediately from the formula that defines it), we get the result. 

Now let < s < t < 1 be general values. Since u is continuous and since 
dyadic rationals are dense, u{t) — u(s) lies in the closure of K, which is equal 
to K U {0}. But since, by Corollary pi u is injective, the difference cannot be 
equal to zero. □ 

We now need to introduce the projective equivalents of the relevant objects. 
Let 7r : H \ {0} — > P 1 (M) be the canonical projection; we let K := ir(K) and 

Mi : P 1 (M)^P 1 (M) 

tt(x) ^ K(Mi(x)) 

(see Figure [7]). Clearly, K is stable by the Mj. 

Proposition 5. Let s £ [0, 1]; let s = 0, a\a2 ■ ■ ■ be its upper binary expansion. 
Then for all u £ K , we have 

lim n(u(t) - u(s)) = lim M ai ... an (u ), (4) 

t— >s n— foo 
t>s 

that is, both limits always exist and are always equal. 
As for lower expansions, the formula becomes: 

lim 7t(u(s) - u(t)) = lim M ai ... a {u ). 

t— >s n— >-oo 
t<s 
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K 



M 00 (K) 
M i(K) 



M (K) 



K 



M W {K) 
M n (k) 



Mi(k) 



Figure 7: Illustration of Proposition |H] 



Proof. The proof is quite similar to the proof of Proposition [2j Just as before, 
the sets M ai ,,, an (K) are all closed and form a decreasing sequence. To ensure 
that their intersection contains a single point, it is enough to check that their 
diameters tend to 0; to check this, it is enough to show that the Mi are con- 
tractive. In the usual metric (where the distance between two points is the 
angle between the corresponding lines), this is not the case: the Lipschitz con- 
stant is equal to 1. But we shall use another metric — the one defined by 

IN 

I 



d(-K(ui),TT(u 2 )) 



bi 



b-. 



where we let Ui = aiV + biW with v 



and 



1 



w 



Informally, what we do is to project on a line rather than on a 



circle: see Figure [7] Since there exists a neighbourhood of the point at infinity 
which does not intersect K, this metric is indeed defined on K and induces the 
right topology on it. Then a few lines of calculation show that in this metric, 
both M and Mi have Lipschitz constant | . Hence we can write 



oo 

e K, { lim M ai ... an (u )\ = f| M ai ... an (K) 

In— >-oo J 1 1 

71=0 



Now let n € N. Then for any t sufficiently close to, but greater than, s, the 
first n bits of the binary expansions of t and s coincide (this is true because 
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we use the upper binary expansion of s), so that we have s — h ai ... an (S) and 
t = h ai ... an (T) for some S,Te [0, 1]. By LemmafJ] we have 

n(u(t)-u(s))=M ai ... an (n(u(T)-u(S))) eM ai ... an (K), 

hence the result. 

As for the symmetric case, the proof is completely analogous. □ 
Corollary 6. The quantity 

Du(s) := lim n(u(t) — u(s)) sgn(i — s) 
t— >s 
*6[0,1] 

is a well-defined, injective, continuous function of s. 

In geometric terms, Du(s) is the direction of the tangent to the curve which 
is the image of u. This result may be interpreted as saying that this curve is C 1 
and convex — which is after all not so surprising if we think of it as the shape 
of the side of an elastic Sierpinski triangle stretched on three nails (see Figure 
[3]) . Note that the parametrization of the curve is much less regular than that 
— as we shall see, it is not even always differentiable. 

Proof. Being given two sequences (<!;);>! and (a^)i>i with values in {0, 1}, we 
claim that 

00 00 

f) M ai ...a n (K) = Ma^jK) 
n=0 n=0 

if and only if both sequences are binary expansions of the same number. The 
proof simply follows, mutatis mutandis, the proof of Corollary [3] we need the 
fact that the Mi are injective and that Mq(K) and M\{K) intersect at a single 
point, which are easily seen to be true. Combined with the previous Proposition, 
the "if part shows that the function is well-defined and the "only if" part shows 
that it is injective. 

As for the continuity, let e > 0. Since Mo and Mi are contractive (in a 
suitable metric), we can find a number n such that the diameter of every possible 
M ai ... an {K) is less than §. Now let s,s' € [0,1] such that \s - s'\ < 2~ n . 
Let s = O.Q4 . . . a n . . . and s' = O.a^ . . . a' n . . . be their binary expansions; let 
w = a\ . . . a n , so that Du(s) £ M W {K), and similarly define w' . Then two cases 
are possible: either w = w' , or w and w' are two consecutive words. However, 
knowing that Mo(ir(wo)) = ir(wo), Mi(-k{wi)) — n(w\) and Mi(tt(wo)) = 
Mo(it(wi)) , a trivial induction shows that in the latter case, M W (K) and M W {K) 
have a common point. Hence the distance between Du{s) and Du(s') is less 
than s, which finishes the proof. □ 

5 Local behavior of u and u a c : rate of change 
5.1 Holder exponent formula 

In the previous section, we have seen that as t tends to s, the direction of 
u(t) — u(s) behaves rather well. We shall now see that its norm behaves much 
more erratically. A naive approach would be to simply try to calculate the 
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derivative u'(s); however, we shall see in Proposition 16 that whenever this 
derivative exists, it is equal to either or oo, so it is not the most relevant 
parameter here. 

We introduce the following parameters, that will provide a better description 
of the local behavior of u: 

Definition 4. Let s e [0, 1]. 

• The upper Holder exponent of u in s is 

ln\\u{t) -u(s)\\ 



a sup (u,s) := hmsup ■ 

t-ys m t 

te[o,i] 



• The lower Holder exponent of u in s is 

ln||u(f) -u(s) 



ctinf (u, s) := liminf 



t-^s In It — s| 
te[o,i] 



• When both are equal, we shall simply call them the Holder exponent of u 
in s, and we shall write 



ln||u(t) -u(s)\\ 

a{u, s) := lim 



t-r+a, In It — s| 
te[o,i] 1 1 

Once we know the Holder exponent of it, we can usually calculate its deriva- 
tive, and then the derivative of u a c b . Indeed: 

Lemma 7. Let s G [0, 1]. Then we have : 

• a;„f(u, s) > 1 w'(s) = ==>■ a in f(u, s) > 1; 

• a sup (u, s) < 1 => u'(s) = ex) =^ a sup (-u, s) < 1, 

where the latter equality is shorthand for "lim f s lffl*2^wll — go 

Proof. This is an easy exercise in real analysis. □ 

Lemma 8. Le£ $ = (a b cj be a linear form, and let s G [0, 1] such that 
Du(s) <£ ker<!>. Then we have: 

. K c h )'(s) = iffu'(s) = 0; 

• (w a c h)'(s) = oo iffu'(s) = oo; 

• (u a c b )'(s) is undefined iffu'(s) is undefined. 

Proof. This is basically an application of the formula for the derivative of a 
composite function. □ 

Note that by Corollary [6] for a given value of this lemma applies for 
all values of s except at most one. In particular, it is easy to check that for 
the functions 4>, ip, x an d £, introduced in pQ, this lemma holds with only two 
exceptions: x'(0) and £'(1) (and in |T|, it is shown that both are equal to 0.) 

Here is the formula that allows one to calculate the Holder exponent: 
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Proposition 9. Let s G [0,1], s — Q.a\CL2 . . . its binary expansion. Then we 
have: 

of m f (u, s) = hmmf — yyv — ; 

n->oo n In ( i ) 

- (5) 
/ \ ,. In M 0l ... r II 

a sup (u, s) = hmsup — 

n-yoo filn 5 



y2) 



a ^ 

/ = I i ) , where the coefficients a, b, c and <i are all nonnegative; 



(compare this with ([3]) and Q). Before proving it, we shall need the following 
lemma: 

Lemma 10. Let f be a linear automorphism of K 2 , let C be a salient closed 
cone that is stable by f, and let C C intC such that C U {0} is also a closed 
cone. Then for all x G C , the ratio 

11/0*011 

II/IIPII 

is positively bounded by constants that depend only on C and C . 

Proof. Let us choose a basis in which C becomes the upper right quadrant. We 
may now write: 

• C = {(x,y) G M 2 | x > 0, y > 0}; 
b 

v c d 

• C = y) G C \ {0} m < | < Af|, where to and M are some finite 
positive bounds. 

Let x = (x, y) G C. Let us use the supremum norm: we have ||a?||oo = max (x, y), 
ll/lloo = max (a, 6, c, d) and ||/(5)||oo = max (ax + by, cx + dy). But we have, 
on the one hand 

max(ax + by, cx + dy) < 2 max (a, b, c, d) max (x, y), 

and on the other hand 

uiax(ax + by, cx + dy) > max (a, 6, c, <i) min (x, y) 

> fimax (a, b, c, d) max (x, y), 

where [i = min (m,M _1 ) depends only on C and C as announced. □ 

Let K' = {awo + bw\ \ a, b > 0} (on Figure ^ it is the cone that lies between 
wq and iZ?i, not represented to avoid overloading the picture). Then it is easy 
to check that it is stable by Mq and Mi, hence also by all of their products. 
By applying the lemma on C = K' and C = K (watch out for the slight 
inconsistency in notations), we get the following estimation: 

Corollary 11. For any finite word w on the alphabet {0, 1} and for all uq G K, 

the ratio 

\\M w (uq)\\ 
l|M»llll(*>)ll 

is positively bounded. 
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Proof of Proposition^ Now let s £ [0,1], and let s' £ [0,1] be a point close 



to, but not equal to, s. Let n 



si < 1. Let 



In \s —s\ 

s = O.ai . . . a„ . . . and s' — O.a^ . . . a' n . . . be the respective binary expansions. 
Let us write w — a\ . . . a n and w' — a' x . . . a„; we then have s — h w (S) and 
s' = h w >(S'), for some S, S' £ [0, 1], hence: 

In ||u(s') — u(s)\\ In u^/i^/ (S')) — u(h w (S)) || 



ln|s'-s| nln(l)+0(l) 

ln\\M w ,(u(S')) -M w (u(S))\ 



(6) 



nln(i)+0(l) 
Now two cases are possible: 

• Suppose w' — w. Then by the previous Corollary, we have 

ln\\M w ,(u(S')) -M w (u(S))\\ =\n\\M w (u(S')-u(S))\\ 

= In HAl^l + In||u(5')-«(5)|| + 0(1). 

(7a) 

• In the other case, by definition of n, w and w' are necessarily binary 
representations of two consecutive numbers. Without loss of generality, 
we may suppose that s < s'. Then we have h w >(0) = h w (l), hence 

\n\\M w ,(u(S')) - M w (u(S))\\ = In 11^(^(50-^(0)) + M w (u(l)-u(S))\ 

Now it is easy to see that, given two vectors u, v that can take values 
in a salient cone, the ratio between \\u + v\\ and + is positively 
bounded. Since both terms in the above sum lie in K, we can write, using 
once again the previous Corollary: 

\n\\M w :(u(S')) - M w (u(S))\\ = 

= ]n(\\M w ,(u{S')-u(a))\\ + \\M w (u(l)-u{S))\\) + 0(1) 
= max(]n\\M w ,(u(S')-u(0))\\, In \\M w (u{l) - u(S)) \\) +0(1) 
= max^ln||Mu,/|| +ln||w(S") -u(0)|[, In +ln|[u(l) - u(S)\\^ + 

Now note that we have 

M w (v ) = Mw'(vi). 

Indeed, we may write w = do • • • OfcOl ... 1 and w' = a® . . . a^lO ... (where 
both last groups of bits have equal, either zero or nonzero, length); the 
identity then follows from the equalities M (v ) — ^v , M 1 (t/ 1 ) = and 
M (vi) = Mi(vq), that are straightforward to check (see also Figure |6" 
Since t/o.i £ K, we get 

ln||Mj|+lnK||+0(l) = ln||M^||+ln|K||+0(l), 
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or in other words, since \\vq\\ = \\vi\\, 

hx\\M w >\\ = ln||M w || + 0(l). 
This allows us to write 

ln\\M w ,(u{S')) - M w (u(S))\\ = (7b) 
= ln\\M w \\ + max(ln||<it(S") -w(0)||, In ||«(1) - + 

Let us now check that in both cases, the middle term is bounded. Since [0, 1] 
is compact and u is continuous, it is clearly bounded from above. For the other 
inequality, note that, in the first case, we have (S" — S) = 2 n (s' — s) > \ > |; 
and in the second case, we have (S' — 0) + (1 — S) = 2™(s' — s) > |, hence 
either (S' - 0) > \ or (1 - S) > j. Now consider the set 

, 1 

l s i ~s 2 \> - 



X:={(si,s 2 )e[0,lY 

and the function 



X — > E>o, 
(si,s 2 ) i ^ ||u(si) - u(s 2 )||. 

Since it is continuous and X is compact, it reaches its minimum value. But by 
Corollary[3] it cannot vanish; hence it is positively bounded. Thus ||u(5")— tt(5')|| 
in the first case, and max(||u(S") — u(0)||, \\u(l) — u(S)\\) in the second case, 
are positively bounded. 

This allows us to merge (l7a| with (|7b|, and to write: 



]n\\M w ,(u(S')) -M w (u(S))\\ = In 11^11+0(1), 
and if we plug this into ([6]), the result follows. □ 

5.2 Case of rational points 

Theorem 12. Let s € [0, 1] (1 Q. Then we know that the binary expansion 
of s = 0,aia2... is eventually periodic; let n be the length of its period, and 
p = Ofe . . .afc + „_! its period. Then M p has distinct, positive, real eigenvalues, 
and we have 

a{u,s) = — 

nm 2 

where X p is the bigger eignevalue of M p . 

In practice, since we know the determinant of M p , we can calculate X p by 
the formula 

X P = \ (x*M p + M p ) 2 - 

Note that Proposition [7] also allows us to infer u'(s). Indeed, we shall see in 
the next section that we never have a(u, s) = 1. 

To prove the theorem, we shall need the following lemma: 
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Lemma 13. Let s € [0, 1], let w = a\ . . . a n be a word on the alphabet {0, 1}. 
Then we have: 

aini(u,h w (s)) =a in f(u, s); 



(u, h w (s)) = a sup (u,s). 



Proof. Let O.6162 • • ■ be the binary expansion of s; let uo e K. It is well-known 
that ^%|r^ is bounded when u varies; since M w is invertible, the reciprocal of 
that quantity is bounded as well; hence the quantity 

ln||M ai ... afcbl ... 6 „('Uo)|| - In ||M bl ... 6n (u )|| 

is bounded. The result follows immediately from Proposition [9j □ 

Proof of Theorem \12\ The lemma we have just shown allows us to drop the 
preperiod of the expansion of s: we may suppose that it is periodic. From the 
formula Q it follows that we have: 

ln||M fe || 

a ini (u,s) = a snp (u,s) = lim 



fcnln (i 



Now since K is stable by M p , we can also say that K is stable by M p . Since 
M p is contractive, this means that it has exactly one fixed point in K, which 
means that M p has a unique eigenvector v p lying in K and associated with a real 
positive eigenvalue A p (we will justify in a minute that it is indeed the bigger 
one). Let /i p be the other eigenvalue; since det(Aip) = (^)™ > 0, it is also real 
positive, and thus the associated eigenvector w p cannot lie in K. Using once 
again the fact that M p (K) is a proper subset of K and considering the positions 
of the relevant vectors relative to each other and to K, it can easily be seen that 
A p > n p . It follows that || ~ A^", hence the result. □ 



Example 1. It follows from Theorem [12| that whenever s is dyadic rational 

In 



In — 

we have a(u, s) = -. — f ~ 0, 737 (indeed, in this case, X p is simply the bigger 



eigenvalue of Mo or Mi, namely |). In particular, we have in this case u'(s) — 
00. 

Note that Lemma [L3| is interesting in itself, even outside the case s £ Q: it 
shows that changing a finite number of bits in the binary expansion of s does 
not change its Holder exponent. For completeness' sake, let us also mention the 
following (fairly obvious) symmetry of u. Let 



P 



(in Figure [5] it simply corresponds to the vertical symmetry axis) . Then 

Proposition 14. For all s G [0, 1], we have u(l — s) = Pu(s). 

Corollary 15. For all s S [0, 1], a- m f (u, 1 — s) = a- ln [ (u, s) and a sup (u, 1 — s) 
OWp(«, s). 
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Proof. Clearly, P 2 = 1, M x = PM^p- 1 = PM P, M = PM X P, and A is 
stable by P. Let s — 0.ai<i2 ... be its binary expansion; then, simply using (j3j, 
we have: 

u(l — s) = u(0.aia,2 ■ ■ .) 

= lim M ai ...a n u 

n— too 

= lim PM ax ... a Pu 

n— >oo 

= Pu(s), 



where a is shorthand for 1 — a and uq 6 A" is an arbitrary starting point. □ 

5.3 Derivative of w 

Theorem 16. For all s € [0, 1], i/(s) is either 0, oo or nof defined at all. 

Proof. Assume the opposite: let s be a number in [0, 1] such that u'(s) exists, 
is finite and positive. Then we have, in particular: 

*->■« \t - s\ 

for some finite positive constant C. 

Let s = 0.aiO2 ... be a binary expansion. Let t n be a sequence, taking values 
in [0,1], such that for all n € N, the binary expansion t n coincides with that 
of s in the first n places. Obviously, such a sequence converges to s. We then 
have, for all n € N: 

\\u(t n ) - u(s)\\ _ ||M 01 ,.. 0n (u(T n ) -u(S„))\\ 
\tn ~ s\ ~ 2-"|T„-S„| ' 

where S n = h~^ an (s) = Q.a n+1 a n+2 . . . and T n = h~*_ an (t n ). Since [0,1] is 
compact, we may extract from the natural numbers a subsequence <f>(n) such 
that SM n ) converges to some value ^oo. 

Suppose that < 1 (by reversing the order of s and t, we can similarly 
treat the case Sec > 0). From now on, we shall require that t n > s. By Corollary 
[TT] we know that 

||-fl4i...q„ || || (u(T n ) ~ u(S n )) || 

||M 0l ... On (u(T n ) - u{S n )) || 
is positively bounded. We also know that 

||M 0l ... On (u(T n ) -u(S n ))\\ 
2 n \T n — S n \ 

tends to C, hence it is positively bounded as well. Suppose for a moment that 
t n is such that T n is identically equal to 1. Then for sufficiently large n, SV n ) 
varies in a closed interval that does not contain 1, and since u is continuous and 
injective, is positively bounded. 
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This means that — 2 _^, ( „) is positively bounded in norm, so that it takes 
values in a compact space. Hence from 4>{n) we can extract yet another subse- 
quence, that we shall call ip{n). such that ^-^("o"' conver g es ; let us call 
its limit. Note that, though we used a particular value of t n to prove that this 
construction was possible, the actual objects that we have constructed do not 
depend on it. 

Now let T be any point in (6*00, 1]; we set t n — h ai ._ an (T), so that for all n, 
T. n = T. We then have: 

|K^(n)) " U( S )\\ _ \\M ai ... a ^ n) (u(T) - ujS^n))) \\ 

™ IV(«)- S I ~ ™ 2-^)\T - S^ {n) \ 

= ||Aloo(M(r) - M(.goo))|| 

|T-Soo| 



which means that the latter ratio is equal to C regardless of the value of T. 

This claim is definitely too strong to be true, and we shall soon prove it is 
absurd. Note that we have 

which tends to as n increases. Hence det ( Mro ) — 0, in other words, has 



1 Oq J U, 111 »_»L1H_,1 V^UIUO, 1 VX oq 

rank 1. The claim can then be reformulated as follows: there exists some linear 
form $ such that the ratio 

$(u{T) - ujSoo)) 
T S 00 

is constant. The interval {S OQl T) contains infinitely many dyadic rational points; 
by Corollary [6j we may find in this interval a dyadic rational value T' such that 
<&(Du(T')) 7^ 0. But Example [l] together with Lemma ^ tell us that in this 
case, $(u(T) — u(Soo)) as a function of T has an infinite derivative in T' . So 
must the whole ratio: contradiction. 

Now let s € [0, 1] flQ, and suppose u'(s) is undefined. Let us apply Theorem 



12 it is only possible if a(u, s) = 1. Let p be the period of the binary expansion 
of s, n the length of the period; then we have A p = (|) n - We know that M p 
has determinant (^) n and that one of its eigenvalues is equal to 1. Hence 
tr M p = l + (i)" + (^) n , and tr25 n M p = 25" + 6" + (f )". On the other hand, 
we know that in the canonical basis, 25M and 25Mi are given by matrices 
with integer coefficients, which means that so is 25 n M p , hence tr25 Tt M p must 
be integer. Contradiction. □ 

Corollary 17. For Lebesgue- almost all s E [0, 1], u'(s) = 0. 

Proof. Let $ be a linear form whose kernel does not meet K; then from Corollary 
[6] it is easy to deduce that the function $o«is monotonic. But it is a well- 
known result (see for example |3_, Theorem 8.19) that a monotonic function is 
differentiable (in its usual sense of "has a finite derivative") almost everywhere. 
Since by Lemma [8j $ o u behaves in the same way as u itself, the result follows. 

□ 
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5.4 Numerical values and estimations of Holder exponent 

Here is a table of all possible periods shorter than or equal to 7 (up to cyclic 
permutation), ordered by Holder exponent: 



s 


Period p 


Length n 


5" tr M p 


Approximate value of a 


I 

3 


01 


2 


7 


1.119 


21 
127 


0010101 


7 


388 


1.096 


11 

63 


001011 


6 


175 


1.086 


5 

31 


00101 


5 


76 


1.085 


1 

5 


0011 


4 


34 


1.078 


19 
127 


0010011 


7 


436 


1.072 


11 


0001011 


7 


472 


1.055 


w 

^127_ 


0001101 


7 


472 


1.055 




001 


3 


16 


1.050 




00011 


5 


88 


1.040 


63 


000101 


6 


211 


1.039 


1 
9 


000111 


6 


223 


1.025 


9 
127 


0001001 


7 


580 


1.012 


5 
127 


0000101 


7 


616 


0.999 


21 


000011 


6 


250 


0.997 


7 
127 


0000111 


7 


628 


0.995 


1 


0001 


4 


43 


0.982 


127 


0000011 


7 


736 


0.962 


1 


00001 


5 


124 


0.936 


-¥- 


000001 


6 


367 


0.903 


127 


0000001 


7 


1096 


0.880 








1 


4 


0.737 



When we look at this table, it becomes apparent that globally, numbers 
whose binary expansions contain long clusters of equal consecutive bits tend to 
have smaller Holder exponents, while those whose binary expansions are more 
"mixed" tend to have bigger Holder exponents. The exponent is lowest when 
the binary expansion ends in an unbroken string of zeroes or ones (that is, for 
dyadic rational numbers — see below) , and seems highest when zeroes and ones 
strictly alternate (that is, for | and its dyadic-rational multiples; the author 
does not know how to prove this). The author was able to prove the following 
statement, that partly confirms this intuition: 

Proposition 18. Let s e [0,1], s = 0.aia 2 . . . its binary expansion. Let us 
define 

1 " 

<2inf, SU p(s) := liminf, sup- ^ |a, - a i+1 1 . 

i=l 

We then have: 

In ^ 

a ini (u,s) < —j+d in {(s); 
In 5 

In ^ 

a sup (u,s) < —j+d sup (s). 
In 2 
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Informally, di n f (s) and d sup (s) give bounds for the asymptotic proportion of 
places where the bit changes value; the latter can also be understood as the 
reciprocal of the asymptotic average length of blocs of consecutive equal bits. 



Corollary 19. If d sup (s) < 1 



In | 



0.263 



then it'(s) = oo. 



lni " - 3.802' 

To prove the proposition, we will need a small technical lemma: 
Lemma 20. For any u G K , b€N and i G {0, 1}, we have 



\M?($\ 



> 



1 



Proof. Let us first express K in the bases that diagonalise M,. A quick calcu- 
lation shows that we have: 















K = 


ja«o + bwo 






lavi + bwi 

















Let i G {0, 1} and u = atii + bwi G K. We then have Mi(u) = |aw,; + hbwi. Now 
let n G N; we have (since Vi and Wi are orthogonal): 



\M?{u)\ 



n , _, I i 2 



(|) n a^+(i)"6w 



bw,; 



.In II „ 1 1 2 



112 

M \2n ||, _ 1 1 2 
.1) \\ bw 42 



\bwi 



> 



(IT 



av, 



i - ii 2 

\a>Vi\\ 2 



\bwi 



, 2n 



ii'S.iii 



= (!)' 

> 



1 ( 

3 \a/ 



l + ±3* 



;kid 2 . 



□ 



Proof of Proposition \l 8\ Let uq G -/?; by Proposition [9] we may then write 

ttinf s) 



. In M 0I ... „ 

lim mi — . 1 . 

«->oo nln(^J 



,. . . In M 0l ... 0n (ito) 
limmi ^jt 

n->oo 77, In ( ~) 



We may now decompose the binary expansion of s into clusters of consecutive 
equal bits. Let n G N; we have: 



Ma 



M ai ... M a 



Ml 1 M 1 ^ ... M l a k <"\-\ JlfcW , 

Ct 1 Q'o <-*fc(ir) — 1 a fcfTi) 
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20 



it 



where the values of alternate between and 1 (so that a\ = a\, a 2 
1 — oi, «3 = Oi, . . . , Qfe(n) — o, n ) and all li are positive. From Lemma 
follows that 

\\M ai ...aAuo)\\ 2 > (f) fe(n) (f) n |K|| 2 

hence (don't forget that In (h) < 0): 

In ||M ai ...Q„(uo)|| 2 k(n) In (|) m ||"o|| 2 
nln(|) " n In (§) n 

When n grows, the third term tends to 0, and the first term's lower limit is by 
definition equal to di n f(x), hence the conclusion. □ 

Finding lower bounds for the Holder exponent seems much harder than find- 
ing upper bounds. Here is a trivial one: we obviously have 

We it, Vt G {0,1}, \\Mi{u)\\ < 

ln(i) 

hence the Holder exponent is always greater or equal to j »fj ps 0.737. However, 

this is useless for finding points where the derivative is zero. The author suspects 
that if the binary expansion of s never contains more than two equal consecutive 
bits, then u'(s) = — but it has yet to be checked. 
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